Theory of Macroscopic Quantum Tunneling in High- Tc c-Axis Josephson Junctions 
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We study macroscopic quantum tunneling (MQT) in c-axis twist Josephson junctions made of 
high- Tc superconductors in order to clarify the influence of the anisotropic order parameter symmetry 
(OPS) on MQT. The dependence of the MQT rate on the twist angle 7 about the c-axis is calculated 
by using the functional integral and the bounce method. Due to the d-wave OPS, the 7 dependence 
of standard deviation of the switching current distribution and the crossover temperature from 
thermal activation to MQT are found to be given by cos 27 and -^cos 27, respectively. We also show 
that a dissipative effect resulting from the nodal quasiparticle excitation on MQT is negligibly small, 
which is consistent with recent MQT experiments using Bi2Sr2CaCu20g+5 intrinsic junctions. These 
results indicate that MQT in c-axis twist junctions becomes a useful experimental tool for testing 
the OPS of high-Tc materials at low temperature, and suggest high potential of such junctions for 
qubit applications. 

PACS numbers: 74.50.+r, 03.65.Yz, 05.30.-d 



I. INTRODUCTION 



The phenomenon of macroscopic quantum tunneling 
(MQT) has attracted much attention of experimentalists 
and theorists for many yearsii^i^ Among several works 
on MQT, Josephson junctions have been intensively stud- 
ied. In the current biased Josephson junctions, the states 
of non-zero supercurrent can move to lower-lying min- 
ima of the potential through the potential barrier by 
MQT. Ivanchenko and Zi'lberman showed the possibility 
of observing MQT in such systems.— As was predicted by 
Caldeira and Leggett, MQT is suppressed by the dissipa- 
tion effect^^iS Later, MQT and the dissipative effects on 
MQT were experimentally observed in s-wave Josephson 
junctionsi^i^i^ii^ 

Recently, the MQT theories for s-wave Josephson 
junction3^i^ d^d^d3,i4 j^^Yg been extended to d-wave sys- 
tems J^ii^iiiii^ii^i^Si^ It was claimed that the influence of 
the quasiparticle excitation on MQT is negligible despite 
of the existence of the line nodes which result from the d- 
wave order parameter symmetry (OPS)^^^ Therefore, 
in c-axis Josephson junctions, e.g., intrinsic junctions^ 
or cross-whisker junctionS j^^i^^'^^ the crossover tempera- 
ture T* from thermal activation (TA) to MQT was pre- 
dicted to be quite high. This can be ascribed not only 
to the weak dissipative naturei^iiiii^ but also the large 
Josephson plasma frequency Up of d-wave junctions, i.e., 
T* cx LUpJ^ On the other hand, in the case of in-plane 
d-wave junctions which are parallel to the Cu02 plane, 
e.g., grain boundary junctions^®"^^ or ramp-edge junc- 
tions^ it was found that the zero energy Andreev bound 
states (ZES)^i^i2^'^ give the Ohmic dissipative effect 
on MQT.i^ii^iiS Later, the above MQT theories have 
been extended to explore macroscopic quantum coher- 
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and propose a d-wave phase qubits 
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The TA-related phenomena in d-wave junctions 
were already observed experimentally by many 
groupsM^^^^^^^^ Until recently, how- 
ever, no experiments of MQT have been reported. First 
successful observations of MQT in d-wave junctions were 
performed by Bauch et al^ and Inomata et al.^ using 
YBCO bi-epitaxial grain boundary junctions j^^'^^i^^'^'^i^^ 
and high-quality Bi2Sr2CaCu208+5 (Bi2212) intrinsic 
junctions, respectively. Subsequently, several groups 
have observed MQTi^i^iSi^i^ and micro-wave as- 
sisted MQTi^i^i^i^ in such systems. They reported 
that T* of c-axis (Bi2212 intrinsic) junctions is high 
(0.5~lK)^£i^i^i^ compared with the high-quahty 
s-wave junction in which T* is at most O.SKjS^ This 
result is consistent with theoretical predictions] ^^d^d^ 

In the previous MQT studies for d-wave junctions, the 
advantage of a large gap value of d-wave superconductors 
has been mainly emphasized. On the other hand, in d- 
wave junctions, there appears a new degree of freedom, 
i.e., a directional dependence of the anisotropic order pa- 
rameter. This directional dependence produces many in- 
triguing phenomenai^i22i'^i^ Additionally, angular de- 
pendence of the Josephson current in d-wave junctions 
explicitly reflects the d-wave OPS,^i^i^ Therefore, how 
MQT depends on the relative angle of the lobe directions 
between two order parameters is an interesting problem. 

In the present paper, we study the dependence of the 
MQT rate (the inverse lifetime of the metastable state) 
on the twist angle 7 for the c-axis twist Josephson junc- 
tion (see Fig. 1) by use of the functional integral and 
the bounce method. We also investigate favorable con- 
ditions for the realization of MQT at high crossover 
temperature T* . Recently, it was found that intrinsic 
junction stacks exhibit a gigantic enhancement of the 
MQT rate,^i^iSi^i^i2a In the present paper, however, 
we treat a single junction in order to clarify the effect 
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FIG. 1: (color online) Schematic picture of the c-axis twist 
Josephson junction. 7 is the twist angle around the c-axis. 



of the anisotropic OPS which is peculiar to high-Tf, su- 
perconductors. Note that such single intrinsic junctions 
made of Bi2212 single crystals have been recently fabri- 
cated by use of the Ar-ion etching method.— 1^ We will 
also present supplementary explanations and results of 
the previous paper a^^i^^'^^ in which the quasiparticle dis- 
sipation effect on MQT is mainly discussed. 

The organization of the present paper is as follows. 
In section II, we derive the effective action for d-wave 
junctions and formulate the theory for the calculation of 
T* . In section III, we present calculated results of T* . 
The advantages of c-axis twist Josephson junctions for 
qubits applications are briefly discussed in section IV. In 
section V, a summary of the results in the present paper 
is given. 



II. FORMULATIONS 

A. Effective action for d-wave junctions 

Let us start from a microscopic model of the d- 
wave superconductor/insulator/d-wave superconductor 
Josephson junction, described by the grand canonical 
Hamiltonian, 



(1) 



where Tii and TL2 are Hamiltonians for the d-wave super- 
conductors 1 and 2: 



^1(2) = II / ^1(2). (^) 



2m 



X 5i(2) - r') V'i(2)a' (r') V'i(2)a W , (2) 

where m is the electron mass, ^ is the chemical poten- 
tial and "0 ('0^) is the fermion field operator. In order to 



obtain the anisotropic order parameter, the anisotropic 
attractive interaction g{r — r') has to be taken into ac- 
count. The third term in Hamiltonian ([1]), i.e., 

nT = Y,ldrdr' [i(r,r')^L(r)V2a(r')+h.c.l (3) 



describes the tunneling of electrons between the two sides 
of the junctions, and 



8C 



(4) 



is the charging Hamiltonian where C is the capacitance 
of the junction and Qi(2) is the operator for the charge 
on the superconductor 1 (2), which can be written as 



Qm^eY^ J drV-Jja)^ (r) iAi(2)^ (r) , 



(5) 



The procedure to derive the effective action is the same 
as s-wave junctions. By using the functional in- 
tegral method,— the ground partition function Z for 

the system can be written as 



Z = 'Dipi'D'ipi'D'ip2'Dip2 exp 



1 



drCir) 



, (6) 



where (3 — l/ksT, ipiip) is the Grassmann field which 
corresponds to the fermionic field operator 0(';/'t), and 
the Lagrangian £ is given by 

a i=l,2'' 

We use the Hubbard-Stratonovich transformation2&i2i 

-kSo'''drJdrdr',^,^{r',T)■4,^{r,T)g{r-r'),|,^{r,T)^|J^{r',T) 



exp 



hp 



dr / drdr' 



l$fr r' t)\ 
g{r-r') 



+ 0T (''^ t) V'i {r\ t) $(r, r', t) 



,(8) 



in order to remove the term ip^ in the Hamiltonian 
Ti-lr). This introduces a complex order parameter field 
^{r,r';T). The resulting action is only quadratic in the 
Grassmann field, so that the functional integral over this 
number can readily be performed explicitly. 

Next, we perform the variable transformation (or 
gauge transformation) "'^^'^'^1^^ 

$(ra,rfc,T) = A(R,r,T)exp[i0(R,r,r)], (9) 
$(r„rf„T) = A(R,r,T)exph^0(R,r,r)]. (10) 

Here, A is a real field, and R = (r^ -I- rf,)/2 and 
r = Va — ri, are center of mass and relative coordinates. 
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respectively. We assume a slow variation of the order pa- duce auxiliary voltage field V which coules to the charge 
rameter in space and time, and hence A (R, r, r) = A (r) operator. Then, the partition function Z becomes 
and (f) (R, r, r) = cj) (R, r) are satisfied. We also intro- 



VAiVA2V(j)iV(l)2VV exp 



h 



drndrh 



I'&il 



.gi {ra - Tb) 



2) 



CV^ 



TrlnG^i 



r 



(11) 



Here, G is a 4 x 4 matrix Green's function 

{ra,Ta,rb,Tb) -i{'^a,rb)S{Ta - Tf,) 



G ^ {ra,Ta,rb,Tb) 

where, 



(ra, Tb) S (tq - n) G2 ^ {ra, Ta,rb, n) 



i{ra,rb) 



t [Va, Tb) e^('?^i(R-^)-'^2(R:T))/2 



(12) 



-t* (ra, rb) e-»('^i(R"^)-'#'2(R-^))/2 



, (13) 



J 



^1(2) 



"^^^^0 + if^ {vi(2) • V) To 
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ywi(2) +' 



h d<f)i 
2~d7 



eVi(2 



fa- Ai 



(2) 



6 {ra - rb) S {Ta - Tb){U) 



In the last equation, we introduced the Pauli matrix 
Ti, the identity matrix tq, the pair potential 



_ Ai(2)(r)e-'^i(2)(R.r) 

^1(2) " i Ai(2)(r)e'*i(2)(i^'-) 



the superfluid velocity i'i(2) = — (fi^/2TO)V(/)i(2), and 
^1(2) = +(~)W2. In the following, we assume that the 
phase varies slowly in space, namely f 1(2) ~ 0. The func- 
tional integrals over the modulus of the order parame- 
ter field A and the voltage field V are performed by the 
saddle-point method, which lead to the gap equation and 
the Josephson equation, respectively] ^^'^'^1^'^ We also as- 
sume that the tunneling matrix element t is finite only in 
the vicinity of the insulating barrier. Then, the partition 
function Z is reduced to a single functional integral over 
the phase difference <f> — <j)i — <j)2. 



with 



Z ^ \ V(j) exp 



(15) 



Next, we expand In G ^ up to the second order in the 
tunneling matrix element t, i.e., 



TrlnG"^ w Trlng"^ - -Tr(gigi) (16) 



-1 ^ [ G^ {ra,Ta,rb,Tb) ^ 

V G2^ {'^a,Ta,rb,Tb) 

i (ra,rf,)(5(Ta - Tf,) 

iHra,rb)5{Ta~'Tb) 



After some calculations, we obtain 

Scff[(f\ 



2^ J V<j>{T)cxp 



where the effective action ScS is given by 
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.hp (J 
Sc'sW\ = J dr— 



C ( h d(t>{T) 



^ \2e dr 
hp i-hp 



^0 



drdr' 



a{T — t') cos ■ 



{r) 



(17) 
,(18) 

(19) 



I3{t — t') cos ■ 



(20) 



The second term in eq. pop describes the dissipation due 
to the quasiparticle tunneling. The third term describes 
the tunneling of Cooper pairs (the Josephson tunneling) . 
The only difference of the effective action between d-wave 
and s-wave junctions is the dependence of A on the rela- 
tive coordinate r, and is included in the two kernels. The 
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dissipation kernel a(r) and the Josephson kernel /3(t) are 
given in terms of the diagonal and ofF-diagonal compo- 
nents of the Matsubara Green functions in Nambu space, 
denoted by G and 



a(r) 



J2 \t{k,k')\^gi {k,T)g2 (fc',-T),(2i) 



k,k' 



/3(r) = --J2\tik,k'fT^{k,T)T^,{k\-T).{22) 



k,k' 



The Green functions are given by 



Q ik,UJn) 
T {k,UJn) 



"(?^n)'+4+A(fe)2' 

ftA(fc) 

(n^„)2+e| + A(fc)2' 



(23) 
(24) 



where = h?k^ /2m — fi and ?kj„ = {2n + \)tt/ P is 
the fermionic Matsubara frequency (n is an integer) . In- 
formation about the anisotropy of the order parameter 
is included in A(fc). In the case of the cuprate high- 
Tc superconductors (the d^2_y2 OPS), A(fe) is given by 
A(fc) = Aocos20 . Here, 9 is an angle between k and 
a-axis. 



B. Effective action for current-biased c-axis twist 
junctions 

Let us turn to the calculation of the effective action for 
the c-axis twist Josephson junctions (see Fig. 1). We de- 
fine 7 as the twist angle about the c-axis (0 < 7 < 7r/4)i2^ 
Such junctions can be fabricated by using the single crys- 
tal whisker of Bi2212,2SZLiSi2a Takano et al. measured 
the 7 dependence of the c-axis Josephson critical cur- 
rent and showed a clear evidence of the dx2_y2 OPS of 
612212,26*78^ 

In what follows, we assume that the tunneling between 
the two superconductors is described in terms of the co- 
herent tunneling, i.e., 



k(fc,fc') 



(25) 



where fe|| is the momentum parallel to the afo-plane. 
For simplicity, we also assume that each superconduc- 
tor consists of single Cu02 layer, Ai(fc) — Aocos26', 
and A2(fc) = Aq cos 2 (6* 7) (see Fig. [1]). Moreover, 
we consider the low temperature limit {ksT <^ Aq). In 
the case where the Josephson junction is biased by an 
externally applied current /cxt, we have to add an addi- 
tional potential contribution linear in (j) }^^^^ At this level 
of approximation, the effective action ScS of the current 



biased c-axis twist Josephson junction has the form 



hp 



dr 



hp 



M f d(l){T) 



u{4>) 



drdr a{T — t ) cos 



(26) 
(r') 



(27) 



where M — C(fi./2e)^ is the mass and ?/((/)) is the tilted 
washboard potential 



C/(0) = -Ej [cos0(r) +r70(r)] 



(28) 



In this equation, 77 is given by 77 = Icxt/Ic, and Ej 
{h/2e) Ic is the Josephson coupling energy, where 



In = 



2e 



(29) 



is the Josephson critical current. In the derivation of the 
Josephson term, we adopted the local approximationJ^^i 
Details of this approximation and its justification are 
given in Appendix \^ 

By substituting Eq. (|22p into Eq. P5)) . we can obtain 
the 7 dependence of Ic, i-C, 



dy 



y/y^ + cos2 26 y/y^ + cos2 2 (6* + 7) 

(30) 



where A'o is the density of states at the Fermi energy. 
Especially at 7 = 0, we get a well-known result^^ 



lc{0) 



2e 



(31) 



Figure [2] shows the twist angle 7 dependence of the criti- 
cal current Iq. It is approximately proportional to cos 27. 
Note that for 7 > 7r/4, the sign of the current changes. 
Thus, in this case, the 7r-junction is formed. These be- 
haviors are attributed to the d-wave OPS of high- Tc su- 
perconductors. 

Next, we will calculate the dissipation kernel air) for 
the c-axis twist junctions. In the case of the c-axis junc- 
tion with 7 = (e.g., single intrinsic junctions),--'- - the 
node-to-node quasiparticle tunneling is always possible. 
Then, the explicit and asymptotic form of a^r) at low 
temperature are obtained as 



a(T) = 



2n^h 



2Tr 



decos^2eKi ^Ao I cos 26* I 



m"^ Rq 1 

lenA^ll'hf 



for 



Ao|r| 



»1, 



(32) 
(33) 



5 



where Ki is the modified Bessel function, Rq = h/Ae^ is 
the resistance quantum and 



(34) 



is the inverse of the normal state resistance R. Equa- 
tion (|33p gives the superohmic dissipation which agrees 
with the previous result&^'^'SS based on the coherent 
tunnehng approximation (Eq. (|25p ). The present model 
shows stronger dissipation effect than the incoherent tun- 
neling model (^(t) ~ and a more realistic tunnel- 
ing model based on the first-principle band-calculation^ 
(a(T) ^ 1/|t|^).^^ Despite this fact, the quasiparticle dis- 
sipation effect on MQT in the present coherent tunneling 
model is still quite weak as will be shown in Sec. III. See 
Appendix B for the derivation of these asymptotic forms 
of the kernel q:(t). 

In the case of nonzero 7, we get 



a(T) 



U|2 /U2 a2 p27T 

' ' ° ° / d0|cos20cos2 (0 + 7)1 





X iCi ( l^Ao|cos26'|j i^i ri^Ao|cos2| 



7)1^ 
(35) 



The asymptotic behavior of ^(t) for nonzero 7 is numer- 
ically estimated as an exponential function, i.e., 



a(r) exp 



■ sin (27) 



Ao|t| 



(36) 



to the suppression of the node-to-node quasiparticle tun- 
neling by the finite value of yL^i^ and is similar to that of 
the conventional s- wave junctions with the constant order 
parameter A: a(T) ~ exp (— 2A|T|/ft) JiiH For nonzero 
7, if the phase (/)(t) varies slowly with time on the scale 
given by H/Aq, we can expand 0(r) — (f>{T') in Eq. ([27l) 




Y/t: 



FIG. 2: The twist angle 7 dependence of the critical current 
Ic- Broken line is proportional to cos 27. 

with respect to r — t' , which results in 



5diss[(/>] 



SM{"f) 



dr . 

V dr 



(37) 



for Ao|T|/?i ^ 1. This exponential behavior is ascribed where 



SM (7) 



drr^a (t) 



fi^\t\^N§ f^'' d9 [cos'^ 26 + cos'^ 2 {d + -f)]E(k)- 2 \cos2d cos 2 {6 + -f)\K{k) 
'0 



A, 







27T 



[|cos26i| + |cos2(6'-f 7)1] [|cos26'| - |cos 2 (6I 7) |]' 



(38) 



for 7 ^ with k = [|cos26'| - 
|cos2(6'-|-7)|]/[|cos26l| -f |cos2(6i4-7)|]. Here, E{k) 
and K{k) are complete elliptic integral of the second 
and the first kinds, respectively. Hence, under the above 
condition, the dissipation action 5diss acts as a kinetic 
term so that the effect of the quasiparticle dissipation 
results in an increase of the mass, i.e., M ^ M + SM 
(the mass renormalization) . This result indicates that 
the influence of the quasiparticle dissipation for nonzero 
7 is quite weaker than that for 7 = 0, and hence it also 
gives negligible contribution to MQT as will be shown 
in the next section. It should be remarked that the 
increase of the mass SM is inversely proportional to Aq. 



C. MQT rate 

We will calculate the 7 dependence of the MQT rate at 
low temperatures. Note that, in the previous studies of 
MQT in d-wave c-axis junctionSfi^'ii'i^ only the limiting 
cases of 7, i.e., 7 = and 7r/8, have been considered. 
The MQT rate is defined by the formulal.>^i2& 

2 

r = lim - ImlnZ. (39) 

/3-+00 f3 

In order to determine F, we employ the bounce approx- 
imationiii^i2& When the barrier is low enough for the 
MQT to occur but still so high that the bounce approx- 
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imation can be valid, T is given by 



(40) 



where Sb = 5eff[</'B] is the bounce exponent, which 
is the value of the the action Scs evaluated along the 
bounce trajectory (I)b{t). Assuming that /oxt is close to 
/c, we can approximate the washboard potential ?7(0) 
as a quadratic-puls-cubic one. Then, we obtain for 
7 = Oi^iiias 



r(0,77) « ro(0,77)exp 
X (l-?7 ) 



-B(0,?7) - Co— ^ 



and for 7 7^ 

r(7,77) K. ro(7,r?)exp[-B(7,r?)] 



In Eqs. gl]) and dUl), 7/ = h^t/Ic{l), Cq 
dxix^/sinh^inx)] ln(l + x'^) w 0.14, 



2e C 
(41) 

(42) 
(277r/8) 



and To (7, 77) is the decay rate without the quasiparticle 
dissipation, i.e.. 



ro(x.) = 12.,(7,.)J^^g[^exp 



36t/o(7,r?) 



5Sa;p(7,77) 



(44) 



where barrier height Uq and the Josephson plasma fre- 
quency ujp are given by 



TT f ^ f^Icil) f, 2N3/2 



(45) 



For 7 = 0, 6M{0) is the mass increment due to the high- 
frequency components {lu > ujp) of the dissipation kernel 
(Eq. 1221)) and is given by 



SM{0) ^ 



TT^Ao 



dx- , — 

VT 



x^ Jo 



Aq 



dss'^Ki{sx) 



(47) 



The second term in the exponent of Eq. (|¥T|) results from 
the low- frequency components {uj < cjp) of the dissipa- 
tion kernel a{T) (Eq. H^ ). 



D. Switching current distribution and crossover 
temperature 

The switching current distribution P(r/) is related to 
the MQT rate r(7, r?) as^^ 



Piv) = -r(7,'?)exp 

V 



1 



(48) 



where v = \dT]/dt\ is the sweep rate of the external bias 
current. The mean value {rj) and the square mean value 
(77^) of the switching current are respectively expressed 
by 



(77) = / dr]'P{f]')f]', 
Jo 

{rf) = / dr,'P{r]')j]''^. 



(49) 
(50) 



Then, the standard deviation (t„ of the switching current 
distribution P{ri) is defined by 



(51) 



In actual MQT experiments, the temperature T depen- 
dence of cr„ is measured. At high temperature, the TA 
decay dominates the escape process. Then, the escape 
rate T is given by the Kramers formula^"-*^ 



r(7, V) = — ^ — exp 



Uo{l,v) 



(52) 



Below T*, the escape process is dominated by MQT. The 
crossover temperature T* is determined from an inter- 
section point of cr„ between T-independent MQT and 
T-dependent TA process. In the low-dissipation (under- 
damping) cases, T* is approximately given b y^^'^°i^^ 



r*(7) 



27rfcB 



(53) 



III. RESULTS 



In this section, we will numerically calculate the MQT 
rate F and the crossover temperature T* for the c- 
axis twist junctions. In the following we choose v — 
42.4mA/s//c(0), /c(0) = 48.54/iA, C = 76.26fF and 
Ao = BOmeV unless explicitly mentioned. These param- 
eters were used in a MQT experiment of Bi2212 intrinsic 
j unctions. 

Figure[3](a) and Figure[3](b) exhibit the 7 dependences 
of the switching current distribution P{i]) and the mass 
increment SM due to the quasiparticle dissipation, re- 
spectively. As 7 increases, P{r]) becomes smeared, and 
the peak position shifts to the smaller values of rj. This is 
because Ic decreases with increasing 7, which makes the 
switching process more likely to occur even for small rj. 
As shown in Fig. [3] (b) , the ratio 6M /M is of the order 
of lO^"' <C 1 for all 7. Therefore, the mass renormaliza- 
tion effect in twist junctions is quite weak, especially for 
large 7. This 7 dependence can be explained as follows. 
As 7 increases, the node-to-node quasiparticle tunnel- 
ing becomes more suppressed. Therefore, the effect of 
the quasiparticle dissipation, namely the mass increment 
SM, is decreasing with increasing 7. 

In Fig. m the standard deviation cr„ of the switching 
current distribution P{ri) is depicted for various 7. The 
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FIG. 3: (color online) The twist angle 7 dependences of (a) 
the switching current distribution P{ni) and (b) the mass in- 
crement 5M due to the quasiparticle dissipation. 
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FIG. 5: The twist angle 7 dependences of (a) the standard 
deviation (7„ of the switching current distribution P{r]) and 
(b) the crossover temperature T* . Broken lines are propor- 
tional to cos 27 and ^/cos2rj in (a) and (b), respectively. 
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FIG. 4: (color online) The standard deviation an of the 
switching current distribution P(r\) for various 7. Doted line 
is proportional to t'^I'^ , Arrows point the crossover tempera- 
ture T* calculated from the approximate formula (|53p . 



crossover temperature T* decreases with increasing 7. In 
the TA regime, the switching current is proportional to 
y2/3 ^87,88 pjg^ [gj .^g show 7 dependences of (a) ct„ and 
(b) T*. They decrease monotonically with the increase 
of 7. Their 7 dependences are approximately given by 
cos 27 and v^cos 27, respectively. The 7 dependence of 
T* can be explained as follows. Due to the d-wave OPS, 
Ic{j) is nearly proportional to cos 27 as shown in Fig. [21 



Thus, we obtain 

T* (xujp(x \/lcil) « \/cos 27 (54) 

if we neglect the 7 dependence of (77). Hence, T* is almost 
proportional to v^cos^. 

For applications to, e.g., d-wave phase qubit,'^^ higher 
T* is desirable for high temperature qubit operation. 
Thus, we will study the dependences of T* on other pa- 
rameters {v, Ic and C) by changing these parameters, 
and clarify the condition for realizing high T* in the case 
of 7 = 0. 

Figure [5] (a) shows the sweep rate v dependence of T*. 
Note that v has to be much smaller than the plasma fre- 
quency ujp in order to apply our theory. The crossover 
temperature T* is an decreasing function of v. How- 
ever, it quite weakly depends on v. More sensitive pa- 
rameters are the critical current /c(0) and capacitance 
C. Figure [S] (b) displays the /c(0) dependence of T*, 
where /c(0) oc Aq (see Eq. ([SI])). As Ic increases, T* 
increases. On the other hand, T* decreases with increas- 
ing C as shown in Fig. [S] (c). The obtained parameter 
dependences of T* are consistent with the simple for- 
mula, T* cx yicTC. In Fig. M (b), we also show the 
result without the quasiparticle dissipation effect. Since 
the ratio SM/M is inversely proportional to Aq, the influ- 
ence of the quasiparticle dissipation becomes stronger for 
smaller /c(0). However, the effect of dissipation is still 
weak even for /c(0) ~ 1 nA, which is consistent with the 
recent experiment with a low Jc Bi2212 surface intrinsic 
Josephson junctioni^ 

Note that, for large Ic, the Josephson penetration 
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FIG. 6: (color online) The dependences of the crossover tem- 
perature T* on (a) the sweep rate v, (b) the critical current 
/c(0) which is proportional to Ao(we set /c(0) = 48.54^A at 
Ao = 30meV), and (c) the capacitance C for c-axis junction 
with 7 = 0. 



Finally, let us comment on the effect of the quasipar- 
ticle dissipation on MQT. As was mentioned above, this 
effect is quite weak even at 7 = 0, in which the node- 
to-node quasiparticle tunneling is possible. In fact, we 
have confirmed that even when the dissipation term is 
neglected, the obtained results are almost the same. Note 
that, in order to see the quasiparticle dissipation effect 
on MQT more clearly, we have to choose at least three 
orders smaller magnitude of the gap Aq as seen from Eq. 
and Fig. [21(b). 



IV. APPLICATION TO QUBITS 

In this section, we briefly discuss the advantage of c- 
axis twist junctions for the qubit application. 

As was shown in the previous section, the quasiparti- 
cle dissipation on MQT is negligibly small. This result 
strongly indicates the high potentiality of c-axis twist 
junctions for d-wave qubit 3*22 193^^95 ^^97^^^^ 

Moreover, in the case of the twist junction with 7 > 
7r/4, the 7r-junction is formed (see Sec. II. B). A supercon- 
ducting ring with the 7r-junction exhibits a spontaneous 
current without an external magnetic field and the corre- 
sponding magnetic flux is half a flux quantum $0 — h/2e 
in the ground state^i^ Thus, the high-Tc superconduc- 
tor ringiS^iiSl including the c-axis twist junction with 
7 > 7r/4 becomes a quiet qubi t^^i^^'^°^'^°^ that can be 
efficiently decoupled from the fluctuation of the external 
field. 

In order to realize the d-wave quantum computer with 
large number of qubits, we have to tune the circuit pa- 
rameter {ujp) of each qubit independently. In actual d- 
wave junctions, however, precise control of LOp, e.g., by 
changing the thickness of the insulating barrier or by the 
oxygen doping, is very difficult. Our results in Sec. Ill 
clearly show that we can artificially and precisely control 
the qubit parameter ujp oc T* only by varying 7 of c-axis 
twist junctions. Therefore, high-controllability of qubit 
parameters is another advantage of c-axis twist junctions 
for realizing scalable quantum computers. 



depth Aj of the junction becomes small. In the case of 
-^c(O) = lO""* A, Xj is typically of the order of 1 /im. In 
order to apply our MQT theory, the size of the junction 
should be smaller than A,/. It should be also noted that 
T* can approximately expressed by 




where Sr and p are the relative permittivity and the 
normal-state resistivity of the junctions, respectively. 
Therefore, T* is almost independent of the size of the 
junction. Thus, in order to obtain high T*, large magni- 
tude of the gap Aq and small magnitudes of the relative 
permittivity Sr and the resistivity p are desirable. 



V. SUMMARY 

In this paper, we have investigated the influence of the 
d-wave OPS on MQT in c-axis twist Josephson junctions 
by use of the functional integral method and the bounce 
approximation. Due to the d-wave OPS, the twist angle 
7 dependences of standard deviation cr„ of the switching 
current distribution and the crossover temperature T* 
are approximately given by cos 27 and y/cos27, respec- 
tively. Therefore, MQT in c-axis twist junction becomes 
a useful experimental tool for testing OPS of high-Tc su- 
perconductors at low temperature regimes. Moreover, 
the influence of the quasiparticle dissipation is found to 
be very weak. This result indicates the high potential of 
c-axis twist junctions for qubit applications. 
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Throughout this paper, we have considered MQT in 
the single junctions where the quasiparticle dissipation 
effect is found to be neghgibly small. On the other hand, 
in multiple stack intrinsic junctions, the long-range ca- 
pacitive coupling between junctionsi22ii2^ gives rise to 
the 0{N^) enhancement of the MQT rate,.^'^'^'^'^'^^ 
where TV is the number of the stack. Therefore, the effect 
of the long-range capacitive coupling on the quasiparticle 
dissipation in multiple stack junctions will be an inter- 
esting subject of future studies. 
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APPENDIX A: NONLOCALITY OF THE (3 TERM 



The effective action in eq. (120)) includes two terms, 
which are nonlocal in the imaginary time. In this paper, 
we have studied dissipation effects caused by the second 
term (the a term) of eq. (PO)) . while the third term (the 
/3 term) has been treated by the local approximation: 
The Josephson coupling is calculated by assuming that 
the local component at r ~ r' is dominant in the double 
integral of the (3 term. In this approximation, we can sep- 
arate the integral by changing variables as T = (r + t')/2 
and t ~ T — t' . Then, the Josephson energy is defined as 

Ej = - dtl3{t). (Al) 



As a result, the Josephson term in the effective action 
takes the usual form 

-Ej / dTcos(/)(T). (A2) 



This approximation is justified in the s-wave supercon- 
ducting junction where the /3-term has an exponential 
form (/3(t) ^ Q-^^o\'^\/''^y In the high-Tc junctions, how- 
ever, the (3 term shows a power-law decay, and the local 
approximation is not justified in general. For example, 
it has been discussed in ref. [s^ that the ohmic power-law 
decay of the /3 term may affect the phase transition of the 
Josephson junctions made of d-wave superconductors. In 
this appendix, we study the nonlocal effect of the /3-term, 
and show that the local approximation is actually justi- 
fied in the calculation of the MQT rate. 

Let us start with the effective action within the local 
approximation, which is given as 



S 



hp 



/a0(7 



V dr 



dr 



^Ec JO 
.hp 

Ej / dT(cos 0(r) -t- 77(/)(t)). 



(A3) 



Here, we neglected the dissipative term. After changing 
the phase variable as ^ = 0— 7r/2+^2(l — ?/), the action 
is written as 



S ^ 



ph 



^EcJo 

+Ej / dr 
Jo 



dr 



d^ 



(A4) 



where Ec = (2e)^/2C. Here, the potential term is ap- 
proximated by a cubic polynomial. We further rescale 
the variables as f = ujqt and ^ = ^2(1 — r])rp, where 
hujo = y/2EcEj{2{l - 77))^/". Then, the action is ob- 
tained as 



h 

^local 



y^(2(l-ry))^/4x5,ocai, (A5) 

(A6) 



df 



1 ~. 



1 7, 



df 
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The same variable change can be done for the general 
action (|^ by using the expansion around the metastable 
state. The result is given by 



nonlocal 



2 V df 



+ df dT'(3{f ~ f) 



6^ 2^ 



where /3(f) = P{f)/Ej and 1/; = (i/;(r) + iIj{t'))/2. The 
normalized action S only depends on the properties of the 
/3-term, and is independent of other parameters such as 
the current, Josephson energy and charging energy. The 
value 5* for the bounce solution determines the prefactor 
in the exponent of the tunneling rate, and is given by 4.8 
in the local approximation. 

In order to deal with the nonlocal /3-term, we need the 
short-time cutoff tq. This cutoff is of order of h/A.^ in 



10 



the d-wave junction considered in this paper, and is much 
shorter than the plasma frequency loq in actual junctions. 
In the following discussion, we focus on the limit ujqTo 
0. For the purpose of this appendix, it is sufficient to 
consider the following simplified kernel (s > 0) (see also 
Appendix B) 



which decays in the power-law form l/|r|'*+^ in the long- 
time limit (r ^ tq), by noting the sum rule / dfp{f) = 1. 

First, we calculate the bounce solution of the action 
including non-local /3 term for the ohmic case (s = 1) 
numerically. We discretize the integral in the action (|A7[) 
as 




-10 -5 5 10 



FIG. 7: Numerically obtained bounce solution for the action 
with the ohmic non-local P term for ti;oro = 1 is shown by the 
solid line. The usual bounce solution for the action with local 
approximation, which is given by ■i/'o('?) =3/ cosh^(f/2), is 
also shown by the broken line. 



E 



1 /i/>. 



At 



At 



At,(A10) 



where At is a time slice and /3ij = j3{Ti — Tj)Ar. The 
stationary condition (55 = is expressed by dS/dtpi — 0, 
and the bounce solution is obtained by the equation 







■ A- 



2Vi 



(Ar)2 



(All) 



Here, we apply the Newton method to solve these non- 
linear equations. For the actual calculation, we intro- 
duce a long-time cutoff Tmax, and discretize the range 
[— Tmax, T"max] with 2N time slices. We consider the 
open boundary condition (V'-jv = V'JV = 0), and choose 
N = 400 and T,„ax = 10 {At = 0.025). 

We show the bounce solution for luoTq = 1 by the 
solid hne in Fig. [T] The calculated bounce solution 
does not deviate so much from the usual bounce solu- 
tion = 3/ cosh^(f/2)) for the local approximation, 
which is drawn by the broken line in this figure. This re- 
sult indicates that the nonlocality of the (3 term is not 
crucial for the shape of the bounce solution. 

We show the action of the bounce solution, 5noniocai 
as a function of the short-time cutoff tq by the solid line 
in Fig. [51 We observe that the value of the action does 
not change so much from 4.8, which corresponds to the 
local approximation. Moreover, in the limit luqTq 0, 
the calculated action approaches 4.8. Similar results are 
obtained also for the general exponent s. 

We can also discuss analytically the irrelevance of the 
nonlocality of the /? term by the perturbative method. 




Perturbation 



0.2 



0.4 0.6 
WqXq 



0.8 



FIG. 8: The normalized bounce action S* as a function of the 
cutoff parameter tq obtained by numerical calculation (solid 
hne) and perturbative method (dotted line). In the limit 
1.^0^0 — * 0, the bounce action becomes 4.8, which corresponds 
to the local approximation. 



We introduce a parameter A, which expresses the nonlo- 
cality of the P term, and write the action as 



^(A) = 5local + AA^ 

AS" = / dfdf'P{f - 
■f ~. 



df 



1 



6 



(A12) 



(A13) 



The local approximation corresponds to A = 0, while the 
action including the full non-local /3 term is obtained for 
A = 1. Here, we develop the perturbation theory with 
respect to A. Up to the first order of A, the action is 
evaluated as 



S ^ SioUMr)] + XASiMr)], 



(A14) 



where V'o('?) = 3/ cosh^(f/2) is the non-perturbed 
bounce solution and 5'iocai[V'o(''')] = 4.8. We show this 



11 



approximate action with A = 1 for the ohmic damping 
by the dotted hne in Fig. [S] As seen in this figure, the 
perturbation method gives a rehable result especially in 
the limit ujqTq — > 0. 

In the following, we show results by the perturbation 
method, which indicates the irrelevance of the nonlocality 
of the P term for a general exponent s by focusing on the 
limit ujqTq 0. We can continue calculation from (jA14p 
as follows, 



For ohmic damping (s = 1), we obtain A5' 
0.084(a;oTo). For s > 2, we have 



AS[iPo{t)] = F 



1 



1 



(A18) 



AS[Mr)] 



1 

16 



df / df'/3(f-f' )(V''-V')^(V'' + V'-2), 



(A15) 

where tp — i^oir) and ip' — V'o ('''')• We can calculate 
the leading correction with respect to small loqTq. For 
< s < 2, we obtain 



AS[Mr)] 



16 



(^0x0)^*7(5), 



(A16) 



where 7(5) is a numerical factor, which only depends on 
s: 



We note that, for any s, the correction by the nonlocal- 
ity of the /3-term disappears in the limit of luoTq 0. 
Therefore, we conclude that the nonlocal effect can be 
neglected at least in the calculation of the MQT rate. 

Finally, we comment on nonlocal effects of the /? term 
in macroscopic quantum coherence. In this case, we have 
to treat the path with many instantons. The interaction 
between instantons governs the coherence of the super- 
position between distinct macroscopic states. In this cal- 
culation, there is a possibility that the nonlocal effects of 
the (3 term become important. Detailed discussion along 
this line is an interesting future problem. 



APPENDIX B: DERIVATION OF THE 
ASYMPTOTIC FORMS OF Q(r) AND /3(r) 



(A17) 



J 



Here, we derive the asymptotic forms of q(t) and /3{t) for c-axis junctions at 7 = q ,3,15^81,85 (jggj-^g order 
parameters as Ai{6) = A2{9) — Aocos20 = A{d). The dissipation kernel q:(t) and the Josephson kernel /3(t) are 
defined as 



k,k' 



(3{t) = --^ |t(fc,fc')|Vi(fc,r).Fl(fc',-r) 
k,k' 

Here, Matsubara Green's functions at the zero temperature are given by 



Gi{k,T) 

g2{k',-T) 

^i(fc,r) 
-Fi(fc',-r) 



e-^^fce(T) 



1 _ a 



El 



El 



^fce(-r). 



e^' fe'e(-T), 



(Bl) 
(B2) 

(B3) 
(B4) 
(B5) 
(B6) 



where Ef^ = y^^^ + A|_, — A (6*) and 8 is the step function. 
Therefore, we have 



fc k' 
+ ^Y.Y.W,k' 



1 



k k' 



El 



E^ 



k' 



E. 



k 



E. 



k' 



1 VVlf '^fc^p-^(^fe+^fcO 



k k' 



(B7) 
(B8) 
(B9) 
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First, we consider the coherent tunneling model where the momentum k dependence of the tunneling matrix element 
is given by tj^ = tSj^ . This approximation is applicable to cross- whisker junctions with a clean insulating barrier. 

Then, the dissipation kernel a(r) and the Josephson kernel /3(r) can be calculated as 



a{T) 



2ui2 r^TT 



Jo Uo 



n 2 



27r2n 

By applying a formula for the modified Bessel functions Ki and Kq, i.e., 

\r\\m\\ 



^0 

poo 1 , , , 



we get Eqs. ([32]) and ([33l>^i^i^ 

a(T) 



and the following expression of /3(t) 



de cos2(26')ii:i 



MAo, 



cos 26*1 



|2 Ar2A2 /•2'r 



Pir) 



27r2?i 



^ d9 cos^{2e)Ko (^Y^" |cos26i| 



|r||A(0)| 



3^2 1 



for T > 



Ao' 



for T>— . 

Ao 



(BIO) 
(Bll) 

(B12) 
(B13) 

(B14) 

(B15) 



Next, we consider the incoherent tunneling model with a constant tunneling matrix: tj^ y = t, which is applicable 
to cross-whisker junctions with an imperfect dirty insulating barrier. After some calculations similar to the coherent 
tunneling model, we obtain^ 



a(T) 



\t\'N§Al 



de |cos 26*1 Ki 1*^°^ 



Rq 1 h 

„ ■ 9 J tor r > -— 



where normal resistance i?i is defined by 



1 47re2 , ,2 , o 



(B16) 



(B17) 



Note that the Josephson kernel (3 (t) disappears due to the angular averaging in the incoherent tunneling model. 

For intrinsic Josephson junctions (c-axis junction with 7 = 0), tunneling matrix element can be modeled as 
tj^ j^/ = t cos^ (29) 5^ y from the first-principle band-calculations, ^4 where the node-to-node quasiparticle timneling 

is inhibited. Within this model, we obtain^ 



and 



air) 



Pir) 



2n^h 



d9 cos'^ (26) K J 



MAo 



cos 2^1 



2TT^h 



decos*{29)Ko 



MAo 



with normal resistance R2 which is defined as 



1 

R^ 



cos 26*1 



\t\'Nl 



45h^ Rq 1 ^ n 

r^ — — r for T> — 

1287rA3 i?2 \rf Aq 



27h^RQ 1 ^ h 

? for T > — 

128^A3i?2 Irl' Ao 



(B18) 



(B19) 



(B20) 



Finally, let us summarize the results of the dissipation kernel a(r) in this Appendix in Table HI including the results 
of s-wave^^ and in-plane d-wave junctionsJ^iSliS^ A remarkable feature of in-plane junctions is the emergence of ZES, 
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which stem from the sign change of the order parameter j ^-'^i^^i'^^'^ In this table, we also show the spectral density 
J{lo) which is defined by the Fourier transformation of a(T) ^^'^°^ 

/>CXD 7 

a{T)oz / ^e-"l"lj(w). (B21) 
Jo 27r 

This table indicates that a wide variety of quantum dissipations can be realized in d-wave junction only by changing 
the junction configuration and the barrier property. 



TABLE I: Asymptotic forms of the dissipation kernel a(T) and the spectral density J(uj) for several types of the Josephson 
junction (JJ). Here, "Exp." denotes a function of the form exp(— a|T|) with a constant a. "Gap type" is defined by its Fourier 
transformation. 



Symmetry 


Model of JJ 




Kernel a(r) 


Spectral density Jiuo) 


Refs. 


s-wave 


JJ without shunt resistance 




Exp. 


Gap type 


13 


s-wave 


JJ with shunt resistance 








13 


d-wave 


c-axis JJ (7 7^ 0) with coherent interlayer 


tunneling 


Exp. 


Gap type 




d-wave 


c-axis JJ (7 = 0) with coherent interlayer 


tunneling 






15^^ 


d-wave 


c-axis JJ (7 = 0) with incoherent interlayer tunneling 


1/r^ 






d-wave 


Intrinsic JJ 






a;* 


85 


d-wave 


In-plane JJ without ZES 








16,97 


d-wave 


In-plane JJ with ZES 




1/r^ 


LO 


16,19^ 


d-wave 


In-plane s-wave/d-wave J J 




Exp. 


Gap type 
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